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Abstract: We show that a recently constructed five-dimensional (5D) model with gauge- 
Higgs unification and explicit Lorentz symmetry breaking in the bulk, provides a natural 
dark matter candidate. This is the lightest Kaluza-Klein particle odd under a certain 
discrete Z2 symmetry, which has been introduced to improve the naturalness of the model, 
and resembles KK-parity but is less constraining. 

The dark matter candidate is the first KK mode of a 5D gauge field and electroweak bounds 
force its mass above the TeV scale. Its pair annihilation rate is too small to guarantee the 
correct relic abundance; however coannihilations with colored particles greatly enhance the 
effective annihilation rate, leading to realistic relic densities. 
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1. Introduction 

The latest years have been marked by the tremendous progresses in observational cosmol- 
ogy. As cornerstones, the detailed maps of cosmic microwave background [|l| and of the 
three-dimensional distribution of galaxies in the Universe |^ have allowed very significant 
improvements in the discrimination among cosmological models and in the determination 
of cosmological parameters. In particular, the case for non-baryonic dark matter (DM) 
as building block of all structures in the Universe has become stronger and stronger: its 
contribution to the present mean energy density is found to be QuMh"^ = 0.105 ± 0.004 |2[ 
(as usual, in this formula the mean DM density is normalized to the critical density 
Pc = 1.879 X IQ-^^h^g/cm^, with h = 0.730 ± 0.019 being the Hubble constant in units of 
lOOkms"^ Mpc~^). The nature of the DM is still unknown. Among viable scenarios, ther- 
mal generation seems the most natural DM production mechanism, and weakly interacting 
massive particles (WIMPs) are among the leading DM candidates: since they are massive, 
their decoupling from thermal equilibrium occurs in the non-relativistic regime; the weak 
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interaction rate with lighter standard model (SM) particles ensures that their thermal relic 
density is naturally of the order of Qbm @1 (for reviews on dark matter candidates, see 
e.g. ii). 

Essentially all theories in extension to the SM predict the existence of new massive 
particles; some of this extra states can indeed be "dark" , i.e. be color and electromagnetic 
neutral, with the weak force (and gravity) as relevant coupling to ordinary matter. A 
better understanding of the mechanism of Electroweak Symmetry Breaking (EWSB) is 
one of the strongest motivations to consider models beyond the SM (BSM); it is indeed 
tempting to search for a framework embedding, at the same time, naturalness for EWSB 
and WIMPs as DM constituents of the Universe. This is plausible whenever there is a 
mechanism preventing the WIMP to decay (or forcing its lifetime to be much longer than 
the present age of the Universe) . The condition of stability is usually fulfilled by introducing 
a new unbroken discrete Z2 symmetry: all SM particles are assumed to be neutral under 
this symmetry, while the WIMP DM candidate is the lightest non-neutral state. Relevant 
examples of BSM theories which aim to resolve or alleviate the SM instability of the EWSB 
and provide DM candidates, include super symmetric and little Higgs theories; in these two 
cases, the Z2 symmetry is identified, respectively, with the R-parity |^] and the T-parity 

Higher dimensional theories may fit as well into this picture: the lighest Kaluza-Klein 
particle (LKP) is potentially a good DM candidate in the class of extra dimension scenarios 
in which a discrete symmetry makes the LKP stable. The simplest models are 5D theories 
with Universal Extra Dimensions (UED) [Q, namely theories where all the SM particles 
are promoted to bulk fields propagating in higher dimensions, where such a symmetry is 
the so-called KK-parity, an unbroken Z2 subgroup of the translation group in the extra 
dimensions Q (see e.g. |jl^ for other frameworks arising from extra dimensions). Despite 
the simplicity of these models, however, UED theories do not shed any light on the EWSB 
mechanism of the SM, whose quantum instability gets actually worse because of the higher 
(cubic) dependence of the Higgs mass on the UV cut-off of the theory. 

One of the main results in this work is to show that stable DM candidates can be 
embedded also in non-universal higher dimensional theories aiming at the stabilization of 
the electroweak scale. For such purpose, we will focus in a recently proposed 5D theory 
in which the Higgs field is the internal component of a gauge field, and Lorentz symmetry 



is broken in the bulk (see e.g. |13| for a brief pedagogical review of such kind of 
models and for further references). Within this framework, a Z2 symmetry (called mirror 
symmetry) has been invoked to improve the naturalness of the model [12|; as a by-product, 
this symmetry guarantees the stability of the lightest Z2 odd particle. Z2 symmetries of this 
kind are less restrictive than KK-parity. Their implementation is particularly intuitive if 
one considers 5D theories on an interval 5VZ2. The mirror symmetry acts on a given field 
and its copy under the symmetry, giving rise to periodic and anti-periodic states along the 
covering circle S^, respectively even and odd under the mirror symmetry. The LKP is then 
identified with the first KK mode of the lightest 5D antiperiodic field in the model, similarly 
to the LKP in UED models, but with the important difference that mirror symmetry is 
not a remnant of a space-time symmetry and hence does not necessarily act on all fields 
in the model. In particular, the mirror symmetry we propose here can be implemented in 
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flat as well as warped spaces, and does not put any constraint on the relation between the 
boundary Lagrangians at the two fixed-points, aside the obvious one of being Z2 even. 

We present here a detailed calculation of the thermal relic density of the LKP in the 
model of |12|. Since Lorentz symmetry in the extra dimension is explicitly broken, there 



is a certain degree of uncertainty in the model mass spectrum. We focus on the region in 
the parameter space where the LKP is the first KK mode of an antiperiodic gauge field, 
roughly aligned along the U{1)y direction in field space. Electroweak bounds require this 
field to be heavier than about 2 TeV, in a range which is significantly more massive than 
the analogous state in the UED scenario as well as most WIMP DM candidates. Since 
the mass is so heavy, the pair annihilation rate for our WIMP candidate is small and would 
tend to lead to the departure from thermal equilibrium at too early times, overproducing 
DM by one order of magnitude or more. On the other hand, the LKP appears within a set 
of other extra antiperiodic fields, most often with the next-to-lightest Kaluza-Klein particle 
(NLKP) being a strongly interacting particle. For reasonable values of parameters in the 
model, the mass splitting between NLKP and LKP turns out to be small, and the NLKP 
becomes the particle triggering the freeze-out and possibly lowering the LKP relic density 
within the observed value (these are known as coannihilation effects In particular, 

the nature of the EWSB in the model implies that typically the lightest Z2-odd fermion is 
the b-, arising from the KK tower associated to the bottom quark. A strongly-interacting 
NLKP gauge boson can be found, instead, in case the mirror symmetry acts on the color 
SU{3)s- For simplicity, we then discuss two classes of viable scenarios: 

1. The LKP coannihilates with the 6„, and gluons are periodic on S^. 

2. Gluons are both periodic and antiperiodic on and the LKP coannihilates also with 
the first KK mode of the antiperiodic gluon. 

Note that in the first scenario there is a further increase in the effective thermally averaged 
annihilation cross section due to a KK- gluon s-channel resonance in 6_ pair annihilations. 
Values of the relic density in agreement with observations are obtained in both scenarios, 
with a moderate degree of fine-tuning (of order few percent), comparable or even lower 
than what one obtains in other cases in the literature when the relic density of the WIMP 
DM candidate is driven by coannihilation effects. 

The structure of the paper is as follows. In Section 2, we introduce the mirror symmetry 



and briefly review the essential ingredients of the model |12], focusing in particular to the 
mass spectrum of the lightest states. In Section 3 we compute the relic density for the 
two scenarios mentioned above and we add some remarks about the fine-tuning needed to 
get the correct relic density. Section 4 concludes. Various details regarding the Feynman 
rules in our model, a one-loop mass splitting computation, the list of all processes relevant 
for the relic density calculation, and the running of the strong coupling constant as are 
contained in the appendices. 

2. Mirror Symmetry and a DM Candidate 

An interesting property of models based on Universal Extra Dimensions (UED) is 
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the possible presence of a Z2 symmetry, remnant of the broken translations along the 
extra dimension, called KK-parity. This symmetry is crucial to make stable the lightest 
KK particle (LKP) and to identify it as a suitable DM candidate [^J. KK-parity inverts 
the segment around its middle point. In terms of a coordinate < y < ttR, it implies 
the invariance of the Lagrangian under the transformation y irR — y. Such invariance 
implies, in particular, the equality of any possible localized Lagrangian terms at y = and 
at y = ttR: Cq = Ct^. Most extra dimensional models which aim to stabilize in one way or 
another the electroweak scale, however, requires Cq 7^ C-,^ and do not respect KK-parity. 



In particular, models based on 5D warped spaces |15| manifestly violate this symmetry. 
It is then desirable to impose some other less constraining symmetry protecting some KK 
modes from decaying. 

The Z2 symmetry we will consider below has been introduced in [^] and allows for 
arbitrary localized terms in the Lagrangian. As it will be clear below, it works for both flat 
and warped spaces. Consider a simple toy model of two interacting 5D real scalar fields 
(jji and (^2 and impose on the system a Z2 symmetry which interchange them: 01^2 ^ 4'2,i- 
Being the Lagrangian invariant under this symmetry, we can impose boundary conditions 
of the following form for 0i and (p2 (in the /Z2 orbifold notation): 

My + '^^R) = My), M-y) = vMy), (2.1) 



where rj = ±. It is convenient to define linear combinations 0± = {cf>i it (/>2)/v 2 which are 
respectively periodic and antiperiodic on the covering circle and with definite orbifold 
parities: <j)±{—y) = ±r](p±{y). Equivalently, one can consider in Eq. ( |2.1] ) the standard 
parity projection (j)i{—y) = rj(j)i{y) instead of 4>i{—y) = r/(/)2(y), resulting in a change of 
parity for (/>_. Under the Z2 symmetry, (j)± — > =b(/)-|-, so we can assign a multiplicative 
charge +1 to 0+ and —1 to (/>_. The localized Lagrangian terms Cq and which can 
include boundary fields as well, can be arbitrary and in general different from each other, 
provided they respect the above Z2 symmetry. We denote such Z2 symmetry as "mirror 
symmetry" in the following. It can also be implemented on gauge fields. For abelian gauge 
groups, it works as before and one is left with two gauge fields, one periodic and one 
anti-periodic. For non-abelian gauge groups, mirror symmetry can be easily implemented 
only when the orbifold twist (or the boundary conditions on the segment) are trivially 
embedded in the gauge group. In such a case, starting from two identical gauge groups 
Gi X G2, the boundary conditions (^]^) leave unbroken in 4D only the diagonal subgroup 
Gj^}- Mirror symmetry changes the sign of all half-integer KK modes, associated to the 
antiperiodic field leaving invariant the integer KK modes of (f)j^. As such, the first 
half-integer n = 1/2 KK mode of (p- cannot decay and is stable. Mirror symmetry acts 
on these fields as KK-parity, provided one rescales R -R/2, but with the important 
difference, as already pointed out, of allowing more freedom in the 5D theory and on the 
localized 4D Lagrangian terms. 

It is straightforward to generalize the action of mirror symmetry for more extra di- 
mensions. For instance, for complex scalars 0i and 02 compactified on a T"^ /Z2 orbifold 

^Notice that the antiperiodic gauge fields A- are not connections of the gauge group G-. The latter is 
not a group, but the symmetric quotient (Gi x G2)/G+. 
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one can have 



(/>i(z + l) = 02(^), + t) = (l)2{z) , M-z)=vMz), (2.2) 

with z properly normahzed dimensionless coordinates on and r its modular parameter. 
As in the 5D case, the lowest KK mode of 0_ is absolutely stable. 

From a model-building point of view, it is of course desirable not to impose mirror 
symmetry ad hoc for the only purpose of getting a stable particle, possibly with the correct 
properties of being a good DM candidate. This is not mandatory but makes the symmetry 
"more natural". In Supersymmetry, for instance, R-symmetry is typically imposed not 
only to have a stable supersymmetric particle but also to avoid baryon-violating operators 
that would lead to a too fast proton decay.^ In the following, in the same spirit, we will 



consider a model [12| where mirror symmetry has been introduced to reduce the fine-tuning 



needed to stabilize the electroweak scale. 



2.1 A Specific Model 



The model we consider is a model of gauge-Higgs unification on a flat 5D space of the 
form R^'^ X /Zi2. It is well known that in models of this sort is hard to get sufficiently 
heavy masses for the Higgs field and the top quark, due to various symmetry constraints, 
including 5D Lorentz symmetry. The latter symmetry, in particular, links gauge and 
Yukawa couplings between each other and does not easily allow to get the correct top 
Yukawa coupling. Due to the radiative origin of the Higgs potential, a large Yukawa 
coupling will also tend to increase the Higgs mass. It has been shown in ^] that by 
explicitly breaking 5D Lorentz symmetry in the bulk (leaving the 4D Lorentz symmetry 
unbroken), one can easily overcome the two above problems of too light Higgs and top 
fields, having now no constraint linking gauge and Yukawa couplings. In the following, 
we review very briefiy the main features of the model — referring the interested reader to 
p^ , p^ , p^ ] for further details — and then consider in some detail the mass spectrum of 
the lightest non-SM states. 

The gauge group is taken to be of the form G x Gi x G2, with a certain number of 
couples of bulk fermions (^'i, ^'i) and (^'2, ^'2), with identical quantum numbers under the 
group G and opposite orbifold parities. We require that the Lagrangian is invariant under 
the mirror symmetry 1^2. The couples (^'i, ^'i) are charged under Gi and neutral under 
G2 and, by mirror symmetry, the same number of couples (^'2, ^'2) are charged under G2 
and neutral under Gi. No bulk field is simultaneously charged under both Gi and G2. 

We can make two different choices for G and Gi^2, depending on whether we double 
the color group or not. We can either take G = SU{3)w x SU{3)s and Gi = U{l)i or 
G = SU{3)w and Gi = SU{3)i^s x U{l)i {i = 1,2).^ As we will see, both choices can 
give rise to a DM candidate with the correct relic density. For definiteness, we focus in the 
following on the case in which Gi = SU{3)i^s x ^^(l)*; the other case can be trivially derived 



^See for a 5D warped model where a Z3 discrete symmetry is imposed to both suppress proton 
decay and have a stable non-SM particle. 

^The doubling of the (7(1) factor is necessary and motivated by naturalness jl2|. 



- 5 - 



in analogy. In total, we introduce (for each SM generation) one pair of couples (^12; ^12) 
in the anti-fundamental representation of /^[/(S)^ and one pair of couples (^121^12) 
the symmetric representation of SU{3)w. Both pairs have C/(l)i,2 charge +1/3 and are in 
the fundamental representation of SU{3)i^2,s- The boundary conditions of these fermions 
and gauge fields follow from Eqs. ( |2.1| ) and the twist matrix introduced in jl^. The 
unbroken gauge group at y = is SU{2) x C/(l) x G+, whereas at y = ttR we have 
SU{2) X U{1) X Gi X G2. We also introduce massless chiral fermions with the SM quantum 
numbers and Z2 charge +1, localized at y = 0. Mirror symmetry and the boundary 
conditions ( |2.1| ) imply that the localized fields can (minimally) couple only to and mix 
with the bulk fermions 

Before EWSB, the massless bosonic 4D fields are the gauge bosons in the adjoint 
of SU(2) X U{1) C 5'C/(3)u,, [/(!)+, gluon gauge fields g+ and a charged scalar doublet 
Higgs field, arising from the internal components of the odd SU{2>)w 5D gauge fields. The 
5C/(3)+_s and SU{2) gauge groups are identified respectively with the SM SU{3)s and 
SU(2)l ones, while the hypercharge U{1)y is the diagonal subgroup of U{1) and [/(!)+. 
The extra U{\)x gauge symmetry surviving the orbifold projection is anomalous and its 
corresponding zero mode gauge boson gets a mass of the order of the cut-off scale A of 
the model (A ~ (3 4)/i2 |l^). The massless fermionic 4D fields, identified with the SM 
fermions, are the zero modes of a coupled bulk-to-boundary fermion system. Differently 
from the bosonic massless fields above, which all have a constant profile along the fifth 
dimension, fermions have a profile which depends on the bulk-to-boundary mixing terms. 
To a reasonable approximation, one can consider all SM fermions localized at y = 0, with 
the exception of the bottom quark, which shows a small wave-function tail away from y = 
and the top quark, which is nearly totally delocalized. All SM fields are even under mirror 
symmetry with the lightest Z2 odd state in the model absolutely stable. Since the bulk 
fermions ^± have 5D Dirac mass terms, in a (large) fraction of the parameter space of 
the model, as we will see below, the lightest Z2 odd state is the first KK mode of the 
antiperiodic U{1)- gauge field, denoted by j4_. 



2.2 Mass Spectrum 

Electroweak constraints fix the compactification scale in the multi-TeV regime. More 
precisely, it has been found in ||l^ that 1/R> 4.7 TeV at 90% C.L. to pass all fiavour and 
CP conserving bounds. The lightest non-SM particles turn out to be in the 1 TeV range 
and thus for all practical purposes we can neglect EWSB effects and consider the mass 
spectrum in the unbroken phase. 

Let us first consider Z2 even gauge bosons. Aside from the massless SM fields consid- 
ered before, we have a standard tower of KK states for all gauge fields, with the exception 
of X, the gauge field of the anomalous U{l)x symmetry, which becomes effectively a field 
with Dirichlet /Neumann boundary conditions at y = O/nR and of Y, the gauge field of 
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the hypercharge C/(l)y, which can mix with X. We have then (n > 1), 

(2n) _ n 

v+ - ^ 



(2"-) 1^ /r, o\ 

= -5 ' (2-3) 



<")=P.^, (2.4) 



where m"^^ and m^^"^ denote the masses of all 5^7(3)^ x U{\)' and 5f/(3)+^s gluon KK 
gauge fields except X and Y . Since Lorentz invariance is broken in the bulk, we have in 
general introduced the Lorentz-violating parameters p and ps, which are the coefficients 



for the gauge kinetic terms of the form for U{1)' and SU{'i)s respectively (see |12] for 



further details). In the following, we will mostly consider the case in which p ~ ps ~ 1, 
the Lorentz-invariant value. When p — 1, the mixing between Y and X is negligible and 
their KK masses are given by 

mi?;) ^ p^ , (2.5) 

(2n) {n - 1/2) 

- P ^ • (2-6) 

The mass spectra of Z2 odd gauge bosons is easily derived, since no anomalies arise here. 
We have 



»?r"=P^^. (2.8) 

The mass spectra for periodic S'[/(2)i-triplet fermions and for all antiperiodic fermions 
is also easily computed, since they cannot mix with boundary fermions. One has 



- ^ J 

where ki are the Lorentz-violating factors entering in the covariant derivative of the fermions 
and Mi are bulk mass terms (notation as in |jl2|). 

The mass spectra for SU {2)l doublet and singlet periodic fermions is more complicated 
and given by the roots of transcendental equations which do not admit simple analytic 
expressions. These equations depend on the bulk-to-boundary mixing terms €\ 2, the 
parameters ki and the bulk mass terms Mi. After EWSB, the SM fermion masses are 
function of these parameters, so that the subspace of the parameter space spanned by 
{ei,ki,Mi) is not totally arbitrary. In addition, the electroweak constraints, perturbativity 
and an estimate of the size of possible Flavour Changing Neutral Currents (FCNC) favours 
a given regime for such parameters. For all quarks and leptons, except the top and bottom 
quarks, el 2 — O-^; — 1- bottom quark we have 2 — 0.2, /cb — 1 and for the top 

quark e| 2 — ^■'^^ — 2.5^. Having fixed e| 2 and ki, the bulk mass terms Mj are derived 
by the known values of the SM fermion masses. 



*This is the only needed and relevant Lorentz violating coupling in the model. 
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Figure 1: Tree-level spectrum for all states with mass < 1/R. The DM candidate is A_' . 

We summarize in Fig. I the masses of the lightest KK states for the typical values of 
the parameters considered above. We report the tree-level mass spectra for both Z2 odd 
and even states for completeness, although the latter do not play an important role in the 
thermal relic density computation. We denote by etc. the first n = 1 KK mode 

of the corresponding antiperiodic fermions ^''L, ^'i and so on. Similarly, for the n = KK 
modes etc. of the SU{2)l triplet fermions. The fields in the fourth and fifth 

column in Fig.|^ (perturbed doublet and perturbed singlet) are the first periodic massive 
resonances of the corresponding SM fields. For ps > p, the lightest Z2 particle is the first 
n = 1 mode of denoted by A^}_\ which will be our DM candidate.^ As can be seen 
from Fig. it does not coincide with the lightest non-SM particle in the model, the latter 
being given by two Z2-even fermions, SU{2)i triplets, which are almost degenerate with 
an other Z2 even fermion, SU{2)l singlet. They all come from the KK tower associated 
to the bottom quark and have a mass ~ l/(5i?). 

Having various free parameters governing the masses of the relevant KK modes, it is 
pointless to compute the mass corrections induced by the EWSB and radiative corrections. 
They can all be encoded in the effective parameters p, ps and ki.^ There is however a case 
in which radiative corrections are relevant and need to be computed. When the n = 1 KK 

^The DM candidate might also be identified with an unstable, but sufiiciently long-hved, particle. In 
extra dimensions, a candidate of this sort might be the radion, whose rehc abundance is typically too large 
p^ . In our scenario, most likely the radion physics will be entangled with the microscopic mechanism 
inducing the 5D Lorentz breaking, which might also provide a stabilization mechanism for the radion. The 
radion physics should then be revised. This analysis is beyond the aims of our paper and may deserve 
further study. 

"As we will see in the following, the region in parameter space where ~ p is the most interesting as 
far as DM is concerned. Strictly speaking, then, we are considering tree-level values of ps and p which differ 
by the correct amount to compensate the splitting induced by quantum corrections. 
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gluons g^J^^ (or KK fermions b^l^) coannihilate with A^^\ the s-channel diagram in which 
a is created in the g^l^ — g^l^ (or h^})—b^^) annihilation might be in resonance and 
amphfy the annihilation in question, decreasing the relic density. Although the absolute 
radiative correction to the mass of (7^^^ or g^^ is irrelevant, being reabsorbed in ps, the 

relative correction matters and it is this the relevant quantity to study — together with 

(2) 

the decay width of — fo'^ quantifying the effect of the resonance. They are also the 
relevant quantities for the if^ annihilation, once the relation between ps and k}j is fixed. 
We have then computed the mass splitting Am^ = 2irL'g}_ — m^^J at one-loop level. Details 
on such a computation can be found in the Appendix B. For the parameter range taken 
above, the result of the splitting is the following: 

Arug = mfl - 2mW ~ _1.4 as mfl , (2.10) 

where ag is the strong coupling constant, evaluated at the energy scale Ps/R- The value 
( |2.10 ) is comparable with the total decay rate Tg of g^\ which at tree-level is purely given 



by the processes 5^^^ qL,RqL,R- For each quark, we get Tg^^/^ = ^ {c-^j^^^^^fasmf}, 
where the couplings c^^/^^^ are given by Eqs.( A.7 ) and ( |A.g| ). Summing over all SM quarks: 



Fg = c^asmf} ~ 1.5 a, mf} , (2.11) 



.g = c asm'g^ ~ i.D as m'g^ 



m 
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where c is the mean value of the couplings c^j^^ squared. As can be seen, T^ ~ |A 
3. Relic Density 

The setup we have introduced is typical for frameworks embedding a cold dark matter 
candidate. There is a tower of massive states which are in thermal equilibrium in the 
early Universe, and a symmetry, the Z2-parity, preventing the lightest of these states to 
decay. We have also shown that it is natural for such stable species to be the A^S^ , i.e. 
a particle which is electric- and color-charge neutral and hence, potentially, a good dark 
matter candidate. As a rule of thumb, the thermal relic density of a massive particle (i.e. 
a particle non-relativistic at freeze-out) scales with the inverse of its pair annihilation rate 
into lighter SM species. In case of the A^^\ we need to take into account that its mass 
splitting with other antiperiodic states can be small: there is a full set of " coannihilating" 
particles, whose density evolution needs to be described simultaneously through a set of 



coupled Boltzmann equations [20, 14, The picture is analogous to what one finds 

for UED models |22, 23|, or sometimes encounters in the supersymmetric frameworks, see 
e.g. |2^, |2^, ^]: if the annihilation rate per degree of freedom of the slightly heavier states is 
larger (smaller) than for the lightest particle, coannihilations tend to delay (anticipate) the 
decoupling of the latter from the thermal bath, and hence to diminish (enhance) the thermal 
relic component. In practice, in the early Universe thermal environment, coannihilating 
states are essentially indistinguishable, since one species is turned into another by inelastic 
scatterings over background particles (these interactions tend to be much faster than pair 
annihilation processes because they are triggered by relativistic background states). The 
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Boltzmann equation governing the freeze out is then conveniently rewritten in terms of the 
total number density n = rij, with the sum running over all coannihilating particles; it 
takes the form pT| : 

— = -3Hn - {(Jesv) (n^ - n^J , (3.1) 

where H is the Hubble parameter, while ngq is the total equilibrium number density, i.e., 
in the Maxwell-Boltzmann regime: 



n^n = —^giui^ K2(^] . (3.2) 



The effective thermally-averaged annihilation cross section {cr^f[v) drives the decoupling 
and reads: 

i^eEv) = Wff ^1 (y) ^^ff ' (3.3) 
with all relevant pair-annihilation channels included in the effective annihilation rate: 




^efT is)-^\l s{s-Aml) ■ ^ ^ ^ 

In the expressions above, Ki{x) is the modified Bessel functions of the second kind of order 
Z; rrii and gi are, respectively, the mass and the number of internal degrees of freedom for 
the coannihilating particle i, with the index i = 1 referring to the lightest state. For all pair 
annihilation processes the kinematics has been written in terms of Peflf ai^d s = A{p'^Q + m\), 
the center-of-mass momentum and energy squared in the annihilation of a pair of lightest 
states; the annihilation process with given initial states i and j needs to be included in the 
effective annihilation rate whenever s > {rrii + mj)^. 



Relic abundances are computed solving numerically the density evolution equation (3.1) 



with the techniques developed in and implemented in the DarkSUSY package |27]. The 
first step is to derive the expression for W^s (s) for generic couplings and mass spectrum 
in the model. For any given set of the free parameters, W^s is then provided as input 
in the DarkSUSY code which makes a tabulation as a function of the effective momentum 
PeS, taking care of resonances and coannihilation thresholds. The Boltzmann equation is 
then integrated numerically in the variable Y = n/s, where s is the Universe entropy den- 
sity; thermal equilibrium Y = Y^q is assumed as boundary condition at the temperature 
T = mi/2, and the evolution is followed up to the point, after freeze-out, when Y settles 
on a constant value: the relic density is obtained just by scaling this value to the value so 
of the entropy density in the Universe today. Contrary to most analyses in the literature, 
we do not perform the computation of the relic density by replacing the thermally averaged 
annihilation cross-section with a truncated expansion in powers of T/mi; our procedure 
gives a more accurate result, especially in case of coannihilation and resonance effects. 

3.1 Minimal DM framework 

We consider first the framework in which the mirror symmetry does not act on the colored 
SU{3)s group, and all gluons are periodic states on S^. In this case, for the typical 
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Lorentz breaking parameter ps and assuming 
as compatification scale \/R = 4.7 TeV. 



set of fermionic parameters introduced in Section 2, the DM candidate A_ shares large 
coannihilation effects with the lightest antiperiodic fermion }P) (see Fig.l); the latter are 
actually two degenerate fermions in the 6 of SU{2>)w (see Table ^ for an account of the 
degrees of freedom of and other relevant particles). The antiperiodic fermions c^^ 
and rl^'* will also be included in the numerical computation of the relic density, although 
their contribution is very small. As in all coannihilation schemes, our results will be most 
sensitive to the relative mass splitting between the DM candidate and the heavier state. 
In what follows we treat the mass of A^}} as a free parameter, or, having fixed the Lorentz 
violating parameter p = 1, use it interchangeably with the compactification scale 
(recall that mA_ = p/{2R)). We also take the mass of b^^^ as a free parameter; this is 
equivalent to introducing a slight departure of the Lorentz breaking parameter kf, from its 
unbroken value kb = 1, having assigned = 0.2, kt = 2.5 and 2 ~ l-^j for all other 
antiperiodic fermions we assume ki = kb and e\2 = 0.1. 

Since electroweak precision tests set a lower bound on the compactification scale at 
about l/R > 4.7 TeV (90% C.L., see [|l2|), the attempt here is to introduce a dark matter 
candidate with a mass of 2.3 TeV or larger. Moreover, A^^^ does not minimally couple to 
the localized fermions, which are the main components of SM fields. The only diagrams 
giving a significant contribution to the A^^^ pair annihilation rate are those with a third 
generation quark in the final state and a third generation antiperiodic fermion in a t- or 
u-channel; this follows from the fact that only for the third generation the antiperiodic 
fermion and gauge boson wavefunctions can have a order one overlap with SM fields. We 
list the set of Feynman rules relevant for this process and the others introduced below 
in Appendix A, while the full list of the diagrams which are needed in the relic density 
calculation is given in Appendix C. In the region of interest for our model, we find that, 
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whenever coannihilations are not effective, tlie tliermal relic abundance of greatly 
exceeds the cosmological limit, see the dotted curve in Fig. ^. 

On the other hand, pair annihilation rates for the b^l^ are much larger and do enter 
critically in the effective thermally averaged cross section: there is a full set of strongly 
interacting final states mediated, in the s-channel, by either the SM gluon or by the first 

(2) 

periodic KK-gluon . In general, it is not relevant to include in our computation states 
with KK number greater than 1; in this case, however, since m^^^ is about twice m^^ 
(recall that m^^^ is of order 1/R, while m[^^ of order 1/(2 i?)), the annihilation diagrams 

(2) 

with the s-channel become resonant and tend to give the dominant contribution 

to the cross section (the effect of resonances induced by second KK particles was first 
pointed out in |25] within the UED context). The enhancement is maximized for splittings 
2m^^ — m'gl which are below the decay width, see Eq. ( ^.11 ), which is indeed much 
smaller than the energy flowing in the s-channel. We find that, on top of the two mass 
parameters and 'm!j^\ the mass of is the third unknown entering critically in our 

analysis; we take it as a free parameter, again in connection to a possible mild variation of 
Pg around its non- violating Lorentz value = 1. 

Finally, there is another relevant issue concerning strongly interacting states we wish 
to mention before going to the illustration of results: we are considering processes taking 
place at a center of mass energy ~ which is about twice the mass of the annihilating 
DM particle. The QCD coupling constant Us should be evaluated at this relatively high 
scale and hence renormalization group effects cannot be neglected, in principle. Indeed, 
the DM abundance is highly sensitive to as which enters quadratically in annihilation 
rates: roughly ^dm oc We have computed the one-loop /3-function for as within our 
framework (see Appendix D for details) and implemented the running numerically in our 
Boltzmann code; at 5 TeV, as turns out to be approximately 0.097.^ 

In Fig. ^ we show the results for the relic density as a function of rrS^ , for a few 

values of the relative mass splitting {m!'^^ — m!~^ )/rrS^ , and taking g^^^ on resonance, i.e. 

rn'gl = 2rn'^\ From the case with zero mass splitting one can read out the cosmological 

upper limit on rrS^ within this framework, namely m'^ < 4.5 TeV, or equivalently the 
bound on the compactification scale 1/R <9 TeV: this scale is comparable to those favoured 
by electroweak precision tests IQ. As expected, the prediction of the relic density scales 
rather rapidly to larger values when the mass splitting among the coannihilating states 
is increased, and, consequently, the value for the mass of the DM candidate approaches 
the region excluded by electroweak physics (in the figure, the light-shaded horizontal band 
correspond to the 3 a preferred region from the combined analysis of data on the CMB 
from WMAP and from the SDSS large scale structure survey [||; models which lay below 
the band correspond to configurations in which ^4^^^ accounts for only a portion of the DM 



^The running of Os was apparently overlooked in Refs. p2| p3| when estimating the effect of coanni- 
hilations of the LKP with strongly interacting KK states within the UED framework. As explained in 
Appendix D, the effect in UED is larger than for our model; since annihilation cross sections scale approx- 
imately as tts X we expect that the values of \/R inferred from the cosmological limit should be 
correspondingly rescaled down. 
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Figure 4: Effective annihilation rate WcS over the center of mass energy squared s, plotted versus 
the effective center of mass momentum p^s- Contributions from single annihilation and coannihi- 
lation channels are displayed. Also shown (dotted line) is the thermal weight function k (units of 
GeV~^ as displayed with the scale on the right-hand side of the plot). 



in the Universe, while those above it are cosmologically excluded). In Fig. ^ we plot the 
thermal relic abundance as a function of the mass splitting of the coannihilating states 
for a model with the minimum allowed compatification radius 1/R = 4.7 TeV, and for 
a few values of ps- In this it can be seen, cosmological constraints restrict the 

Lorentz breaking parameter of SU{3)s roughly in the range [0.9, 1.2]; the interval is not 
symmetric around ps = 1 since, in the Boltzmann equation, annihilations take place at 
a finite temperature. For ps < 1, temperature corrections drive the process at energies 
always slightly above the resonance. In the opposite regime the resonance is always met, 
provided one considers sufficiently high temperatures; if ps is large, the temperature at 
which the resonance is hit is too large compared to the freeze-out temperature and the 
models becomes cosmologically excluded. Curves for the four sample values of ps overlap 
at a mass splitting of about 30%, beyond which coannihilation effects induce negligible 
changes on the relic abundance. 

Conservatively, we include in the relic density calculation all states with a mass split- 
ting below 50%. In Fig. ^ we illustrate better the role of coannihilations and of the Boltz- 
mann suppression when mass splittings become too large. We consider the model with 
1/R = 4.7 TeV, Ps = 1 and ki = 1, with relic density of about 0.1, and plot the effective 
annihilation rate W^s, as defined in Eq. (3.4), over the center of mass energy squared s, 
as a function of the effective center of mass momentum peff- Contributions to VFeff from 
the individual annihilation and coannihilation processes are shown; coannihilations appear 
here as thresholds at y/s equal to the sum of the masses of the coannihilating particles, 
with the NLKP entering first, and CL^ at a slightly larger peff- The threshold effects are so 
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mass splitting and compactification radius in 
models with matching the best fit value 
from cosmological observations, for a few val- 
ues of ps ■ 




Figure 6: For a few selected values of the 
compactification scale, regions in the param- 
eter space ps~niass splitting in which Q.h'^ is 
lower or equal to the best fit value from cos- 
mological observations. 



sharp since coannihilation rates are large, but also because the number of internal degrees 
of freedom for the antiperiodic fermions is much larger than that for (see Table |3|). 
Also shown in Fig. ^ is the weight function K{pef[,T) defined implicitly by rewriting the 
thermally averaged cross section in Eq. ( |3.3D as 

WeSV) = dpeS K{peS,T) . (3.5) 

Jo s 

The function k contains the Boltzmann factors (hence it is exponentially suppressed at 
large Pefr) and a phase-space integrand term (hence it goes to zero in the PefT ~^ limit). 
It can be view as a weight function, since at any given temperature T, it selects the range 
of peS which are relevant for the thermal average. In Fig. ^, k is plotted as a function 
of peS at the freeze-out temperature (defined as the temperature at which Y is equal to 
twice the final asymptotic value) in units of GeV~^ and with respect to the scale shown 
on the right-hand side of the plot. On the top of the panel, the tick mark with the label 
'0' corresponds to the momentum at which k has its maximum, while the tick mark with 
label —n indicates the momentum at which k is 10""" of its peak value. Coannihilation 
effects are relevant if they provide a significant enhancement in the effective annihilation 
rate within the range of momenta in which k is not too small compared to its peak value; 
this is clearly the case for the b^S^ in the example displayed. Also notice that the effect 
induced by the c^"* is not negligible by itself, however it gets marginal when superimposed 
to the one from the b^l\ 

Figs. ^ and ^ summarize the picture within our minimal DM framework. We select 
models whose thermal relic density matches the best fit value from cosmological observa- 
tions Q.DMh'^ = 0.105. As already explained, there are three relevant mass parameters 
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in the model: or equivalently m[^^ or equivalently the relative mass splitting 

between b^}) and , and vri~gl or equivalently ps ■ In Fig. ^ we select a few values of ps 
and find the isolevel curves for O/i^ in the plane of the other two; in Fig. |6|, instead, a 
few values of the compatification scale are considered and correlations between the other 
two parameters derived. The general trends we see here are essentially along the lines we 
have already discussed for Figs. ^ and |3|; we display more clearly the strict upper limits on 
\/R (about 9 TeV), and find that the NLKP-LKP mass-splitting needs to be at about the 
7% level or smaller. The required range of ps, for a given compactification scale and mass 
splitting, is also displayed. 

We have definitely found a tight interplay among the parameters in the model; the 
procedure of embedding a DM candidate in this minimal scenario has been successful, 
pointing to a limited set of patterns in the parameter space. 

3.2 DM in the framework with a copy of SU{?>)s 

If the Z2 mirror symmetry acts on the colored SU{3)s group, the first antiperiodic KK 
gluon g^^\ which has a mass of order 1/(2 i?), enters critically in the computation of the 
relic abundance for the . In most extensions to the SM, strongly interacting gauge 
bosons are the particles with largest pair annihilation cross section per internal degree 
of freedom, hence tend to give the largest possible coannihilation effects. This has been 
verified also in the extra-dimension context studying the coannihilation of the LKP with 



the first KK excitation of the gluon in UED |22 



We discuss the phenomenology in our model referring again to the three mass param- 
eters introduced above. Note, however, that in this case we select values of ps to fix both 
the mass of (we implement the tree-level relation m^g^} = ps/{2R)) and the mass of 
Qj^ (through the 1-loop mass splitting as found in Eq. ( 2.10D ). We start by examining the 



effect of g_ coannihilations alone. In Fig. ^ we set k^, = hi = 1.5, removing all antiperi- 
odic fermions from the coannihilation list, and discuss the effect of degeneracies in mass 
between g^^^ and A^}} . In the limit of zero mass splitting we find as upper bound on the 
compactification scale 1/R < 11 TeV. As expected, the bound on 1/i? found within the 
minimal scenario has been relaxed. We also find that, at the lowest allowed value for 1/R, 
{m'gl — rrS^ )/rn'^ < 6% must old. Even in the present framework, g^^^ coannihilations 
appear as sharp thresholds in the invariant rate. The channels contributing to the anni- 
hilation rate are listed in Appendix C. They include the case of annihilation into quarks 
(2) 

with the i'^ ^ s-channel; however, this process always takes place slightly off-resonance, 
since |Amg| is of the same size as F^, and it is then always subdominant with respect to 
the process with gluon final states (one may compare the behavior of the g^^^-g^^ term in 
the right panel of Fig. |^ as a function of peff, with the analogous for the b'-^-b'-y ter m m 
Fig. ^, where the enhancement due to the resonance is instead evident). 

The general framework, with both g^^ and 6*^^ playing a role in the relic density 
calculation, is illustrated in Figs. ^ and The picture is not a mere overlap of two 
distinct coannihilation effects. As we have already mentioned, at a given 1/R and given 
mass splitting between and , the mass splitting between g^^ and A^S^ sets also 
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Figure 7: Le/t Panel: Relic density versus the A^'^^ mass, for a few value of the relative mass 
splitting between g^Jl"^ and aL^"* . Antiperiodic fermions have been decoupled assuming they have 
a mass splitting larger than 50%. Right Panel: Effective annihilation rate WcS over the center 
of mass energy squared s, as in Fig. |4[ but for a model with yi^^ coannihilations decreasing the 
relic density of ^L^'' to the level of the best fit from cosmological observations. The thermal weight 
function k is shown as a dotted curve; see Fig. ^, and the relative discussion in the text, for further 
details. 



rrig^ and hence whether the b_ pair annihilation is resonantly enhanced or not. The 
second effect is due to the fact that we are superimposing coannihilations from states 
with different annihilation strength, and, especially, different number of internal degrees of 
freedom {g = 24 for g^^): for equal mass splitting, the matching needs to be done at the 
level of annihilation rates per degree of freedom. The net effect can be both of increasing 
or lowering the thermal relic abundance for A^^\ E.g., if we add, on top of a configuration 
with efficient g^^^ coannihilations, a ft^'* state with small mass splitting with the but 

(2) 

with mass significantly displaced from the resonance, we are effectively including a 
set of passive degrees of freedom: maintaining the tower of states in thermal equilibrium 
becomes more energetically expensive, the freeze-out is anticipated and the thermal relic 
density increased. This is what happens at small m^g} — and small — in the 
throat region of Fig. ^ 

Introducing the g^^^ in the framework has enlarged the regions in the parameter space 
which are compatible with the cosmological constraints; still, the tight correlation patterns 
among parameters in the model, which had emerged in the minimal scheme, are maintained 
here, although in slightly different forms. 

This feature might be view as a sign of fine tuning. To better quantify this point, in 
analogy to the study of naturalness of radiative electroweak symmetry breaking |^^, we 
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introduce the fine-tuning measure |3C, 31 1: 



max 



Oln(a) 



(3.6) 



where a labels any of the free parameters in our model. In the minimal DM framework 
of Section 3.1, changes from about 35 in the lower part of Fig. ^ to about 8 for the 
models with largest ps- In the model with antiperiodic gluons, the parameter space with 



small m 



(1) 



m 



(1) 



and intermediate m 



(1) 



- nS^ in Fig. ^ has a minimum A*^ of about 

fine-tuning correlated to the A-'^^^-g^'^^^ mass 
splitting of about 34; finally in the throat region, in which both mass splittings are small, 
the interplay among the parameters reaches its maximum and, correspondingly, A^ can be 
as large as 50. Such moderate degree of fine-tuning (A^ < 10 is expected in a "natural" 
model) is comparable or even smaller than what one obtains in other cases in the literature 
when the relic density of the WIMP DM candidate is driven by coannihilation effects, see, 

e.g., m. 



7, while in the limit of large m 



(1) 



m 



(1) 



4. Discussion and conclusion 



We have shown how to embed WIMP dark matter candidates into non-universal flat higher 
dimensional theories aiming at the stabilization of the electroweak scale. We have focussed 
on the model of [12| and shown that, in a large fraction of the parameter space, the lightest 
antiperiodic particle is a stable gauge field with the correct properties for being identified 
with the DM in the Universe. Although electroweak bounds force its mass to be larger than 
about 2.3 TeV, and its interaction rate is rather small, coannihilation and resonance effects 
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involving colored particles can delay its decoupling from thermal equilibrium, and allow its 
relic abundance to be within the range currently favored by cosmological observations. 

The picture we have introduced is rather unusual, since the WIMP dark matter candi- 
date is significantly more massive than in standard (e.g. SUSY) scenarios, and its coupling 
to the SM is essentially limited to third generation quarks. The phenomenology of DM 
searches for this model is less appealing than in other frameworks; in particular its scat- 
tering rate on ordinary matter is suppressed and mediated mainly by radiative effects 
involving virtual bottom and top quarks. Moreover, its zero temperature pair annihila- 
tion rate (again mainly into bottom and top quarks) is small, at the level of few times 
10^^^ cm^ s^^ (see Figs. ^ and |^, making it hard to detect annihilation signals in dark 
matter halos. 

On the other hand, embedding the dark matter candidate in the model of [^] intro- 
duces favored patterns in the parameter space; tests of this framework at future colliders 
may indeed give crucial information on the DM scenario proposed in this paper. 

It would be interesting to implement the mirror symmetry used in this paper to achieve 
a DM candidate in warped models as well, mainly for the case of gauge and matter fields 
in the bulk. One could for instance double the U{1)y gauge field and identify the DM 
candidate as the lightest mode of A^. Since couplings to fermions might be larger and 
the reference mass scale lower with respect to the flat case, it might be feasible to match 
the correct relic density without coannihilation effects with colored particles. On the other 
hand, in order to make such construction as natural as in our case, it would be desirable 
to find some other motivation to introduce the mirror symmetry. 
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A. Feynman Rules 

In this appendix we give some details about the Feynman rules of our model, focusing in 
particular on vertices relevant for the calculations of Section 3. The Lagrangian (aside 
from ghosts and gauge-fixing terms) is given in Eqs.(2.4)-(2.7) of |12|. The gauge-fixing 
terms (and the corresponding ghost terms) we use are of the form 

-^5/ = -4 E (^M^'^'" -^pd,Atf , (A.l) 

^ a 

for all gauge groups. All cross-sections have been evaluated in the ^ = 1 gauge. Since 
ghosts and gauge bosons A^, A^ are purely bulk fields, ghost, 3-bosons and 4-bosons 
vertices are easily derived from the usual standard ones. One has only to take into account 
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the Lorentz violation in the fifth dimension replacing — > pA^, — > pd^ and take the 
linear combinations (l)± = {4>i ± (f)2)/V^ for U{l)i and SU{3)i^s gauge and ghost fields. 

Fermion-gauge couplings are more involved, due to the non-trivial profile of fermions 
in the extra dimensions. The interactions between a gauge boson KK mode p with fermion 
KK modes m and n can be written as iT'^g4_^^{c\^^'^^ Pl + '^^R^a^'^^ Pr) ■ The coupling (74 
is the 4D gauge coupling, related to the 5D one as (74 = g^/^/TTrR, the indices p,m,n 
run over even (odd) integers for Z2 even (odd) fields and c^]^^^ are the integrals of 
the wavefunctions along the 5th dimension involving respectively left and right fermion 
components and broken or unbroken gauge field components A^a- In terms of the mode 



expansion (see Appendix of |12] for further details) 



'^L/R=Y.^L/R(y^^L/R' 



in) 



(A.2) 



E(n) (n) 
9l/rXl/ 



R ' 



where n in Eg. ( A.2 ) is even (odd) for periodic (antiperiodic) fermions, one gets 



^(m,n,p) 
^L/R,a 



2-kR 



V2^i dyfliiy) /^/^(y)/^^/^) + + ^{"^^^^^(y) 



:,L/R^ 



L/R"^ 



:,L/R'^ 



(A.3) 

where f^l{y) is the wave-function of the p*'* KK mode of A^j^^aiv)- 

As one can check from the Feynman diagrams listed in Appendix C, the relevant 
couplings in our calculation are: 

• p = 0,m = n: only gauge bosons of the unbroken SM SU{2)l x U{1)y x SU{?>)s 
gauge group have zero modes, with a constant wavefunction: = The integral 



in square brackets in Eq. ( A.3 ) is normalized to be 1 in order to have canonical fermion 
kinetic terms: 



^(0,n,n) 



1, 



^L/R,a - (A-4) 

implying universal couplings for all fermions, as expected from the unbroken gauge sym- 
metry. 

• j) = m=l,n = 0: one gets 



.(1,1,0) 
^R,a 

J1,1,0) 



± 



k{k^ MR) e 



V2tt KM {k^ + M^R'^) ^/Z^ 

kj (kj T MjR) ej 
V2ttM,R {kf + M^R^) 



(A.5) 
(A.6) 



In Eqs.(A.5) and ( |A.6| ), the two different signs refers to the two towers of antiperiodic 
fermions with same mass and quantum numbers and the Z factors are those appearing in 
Eq.(2.18) of 1 17 1 taken at a = (no EWSB). These factors are typically ~ 1, aside from the 
top quark where they can be substantially bigger (~ 4 in the chosen setup). In Eq. ( [A.(]| ), 
i = u,d, depending on the doublet component, and M in Eq. ( A.5[ ) should be identified 
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with Mu or M^, depending on the singlet field under consideration. Similarly for e and k. 
Antiperiodic fermion and gauge boson wavefunctions vanish at y = and thus the overlap 
with the SM n = fields is small, 0(e), except for the top and the left-handed bottom 
quark, for which one has an overlap ~ 0(1). 

• p = 2,m = n = 0: we are interested only to the fermion gauge couplings to g^'^\ the 
first KK mode of SU{3)s- One gets 



.(2,0,0) 
.(2,0,0) 



V2 
V2 



1 + 



MR 



1 + 



2 MiR 



coth(^) 



i=u,d 



TTkiikf + AMfR^) 



coth 



IT MiR 



(A.7) 
(A.8) 



This is a KK-number violating coupling, due to the localized Lagrangian term. As can be 

), this coupling is ~ \/2 for all SM fermions, but the top and 



seen from Eqs.(A.7) and 
the left-handed bottom for which it is much smaller (~ \pllZ\ ). 

• p = 2,m = n = l: again, the only coupling relevant for us is the one with g^'^^ . Only 
bulk fields are involved and the computation is trivial, giving 



.(2,1,1) 



1 



(A.9) 



All effects involving KK states with p > 2, with the exception of the possible gluon 
resonance state for p = 2, have been neglected. 

Analogous considerations can be done for the couplings between fermions and the 
would-be Goldstone bosons A^. The vertices can be written as —kT°-g/^^^{d^]^^'^^PL + 



d 



{m,n,p) 
R,a 



Pji) where k is the Lorentz breaking factor and 



d 



(m,n,p) 
L/R,a 



\/2ttR 



2-nR 



p(") 



e(m) 



{n) 



. ftm) 
'2,R/L^ 



(A.IO) 



The only coupling relevant for us is the one with the colored would-be Goldstone bosons 
p = m = 1, n = for which one has 



^(1,1,0) 



m 



(1) 

/- 



.(1,1,0) 
'■L/R,g 



(A.ll) 



B. One-loop Gluon Mass Correction 

One-loop computations on orbifolds are conveniently performed by using the method of 
images to map the propagators on 3^/7^2 to those on the covering circle [32|. In this 
way, the vertices conserve the KK number and the KK violation induced by the boundaries 
is all encoded in a term in the propagator of the bulk fields. 

As discussed in the main text, the only radiative correction of interest to us is the 
mass splitting Amg 



%(2) 



2m^(i). 



There are three classes of radiative corrections: i) 



bulk (finite) corrections induced by bulk fields, ii) localized (divergent) corrections induced 
by bulk fields and iii) localized (divergent) corrections induced by boundary fermion fields. 
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The corrections i) and ii) are one-to-one, in the formalism of |32|, to loop corrections with 
respectively an even or odd number of insertions of KK- violating propagator terms. 

This picture is valid in the limit of vanishing bulk-to-boundary mixing terms (e — > 0), 
that is a very good approximation for all the fermions but the top. In the latter case, the 
calculations are more involved, since ej ~ 0(1) and the corrections ii) and iii) cannot be 
separated. We have nevertheless checked that the effect of e is negligible in the radiative 
correction. Indeed, by taking the opposite limit e — > oo, in which several simplifications 
occur, the top contribution to the mass splitting varies ~ 1% with respect to the e = 
contribution. For all practical purposes, we can thus safely take e = for all SM fields and 
consider separately contributions ii) and iii). 



B.l Bulk Contributions 

Bulk contributions are easily computed. Since there are no bulk fields charged under both 
SU{3)i^s and S'C/(3)2,s, mirror symmetry constrains the one-loop mass corrections to the 
gluons gi and g2 (and hence g+ and g-) to be equal. Divergencies appear but they are 
associated with the renormalization of the 5D coupling constant and the Lorentz violating 
parameter ps. The former does not alter the mass spectrum and the latter dependence 
clearly cancels in computing Anig. What is left is a finite universal correction, similarly to 
the case of ||3^ . The purely bosonic and ghost contributions are as in [^] , once one rescales 
— > Ps/R, since the Feynman rule for periodic and antiperiodic fields are essentially 
the same. Antiperiodic odd fields running in the loop give only rise to a phase (— )"" 
after a Poisson resummation on the KK modes is peformed. The contributions of virtual 
odd fields in the diagrams is proportional to J2w=i(~)^ / '^'^ ~ ~3C(3)/4, and equals then 
(—3/4) times the ones of the corresponding even fields, giving a partial cancellation. In 
total, the gluon and ghost contributions equal 



9 a,C(3) pI 



g.+gh. 8 



4 



(B.l) 



Eq. (B.l) is valid for all periodic (even n) and antipeioridc (odd n) modes and is inde- 
pendent of the KK number of the external gluons, with the only exception of the n = 
massless QCD gluons for which one clearly has <5m-^(o) = by gauge invariance. 

Fermion loops are similarly treated, although now the Lorentz breaking cannot be 
simply rescaled away. For a couple of fermion pairs (^'1,2^1,2) in the fundamental repre- 
sentation of SU{3)i^2,s with bulk mass Mi and Lorentz breaking parameter ki, one finds 



,(n) 



fer. 



Ogk 
^2 



TT' 



E 

w=l 



-2wXi/ki 



1 + i-Y 



l + 2w- 



ki 
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where Aj = irMiR and we have neglected negligible corrections 0(1 — k"^ / pi) in Eq. ( |B.2| ). 
The terms proportional to 1 and (— )'" in Eq. (|B.2| ) correspond (for Z2 even gluons) to pe- 
riodic and antiperiodic fermion contributions respectively. As above, a partial cancellation 
of the mass correction is induced by antiperiodic fields. Again, the mass correction (|B.2| ) 
is valid for any KK number of the external gluons, but the n = gluons. 
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B.2 Localized Contributions from Bulk Fields 



Due to the presence of one non-diagonal propagator, no sum over KK modes has to be per- 
formed in the Feynman diagram loop associated to these contributions. The diagrams are 
effectively four dimensional and logarithmically divergent. Such divergencies are cancelled 
by introducing boundary kinetic counterterms for the gluons at the orbifold fixed points. 
Strictly speaking, this kind of contributions would then be uncalculable, depending on the 
arbitrary renormalization prescription chosen to cancel these divergencies. It is however 
possible to estimate their effect by assuming that they are dominated by the calculable 
radiative corrections of the model. In other words, we require as renormalization prescrip- 
tion the vanishing of these counterterms at a scale of energy equal to the cut-off A of the 
theory. 

The mass correction is encoded in the r]^i, coefficient 11 of the gluon vacuum polarization 
term, taken at = ""^^(,1) • Contrary to the bulk terms, boundary corrections also induce 
mixing between the KK modes, so that a diagonalization of an infinite mass matrix should 
be performed in order to get the mass eigenvalues. All off-diagonal components are however 
one-loop induced, so that at one-loop level we can safely neglect such terms and focus only 
on the diagonal two-point amplitudes. Since the 11 factor is given by a 4D loop diagram, 
its form is the same for periodic and antiperiodic gluons. The only non-trivial issue is 
the sign of the mass correction. The latter is fixed by the boundary conditions ( |2.l| ).^ 
The ending result is that no localized mass term is induced at y = 0, whereas at y = irR 
the periodic and antiperiodic contributions are equal. The localized mass contributions 
induced by gluon and ghost fields is found to be (n > 0) 

where mg{n) = Psn/{2R) is the tree-level mass for periodic and antiperiodic gluons. The 
localized contributions induced by bulk fermion fields vanish trivially because the KK- 
violating terms in the fermion propagator contains a 7^ factor which results in a vanishing 
trace over the spinor indices. From Eq. ( |B.3| ) we get the following one-loop contribution 
to Am^: 



independently of the cut-off A. 



B.3 Localized Contributions from Boundary Fields 

The contributions from colored fermions localized at y = is straightforward. Being a 
purely 4D contribution, it is logarithmically divergent and will be renormalized as described 
before, requiring the vanishing of the localized operator at the scale A.^. Boundary fermions 



^Instead of considering periodic and antiperiodic fields, as usual, one could alternatively consider an 
/{7i2 X Z2) orbifold where all fields are periodic but with different orbifold parities at j/ = and y = tvR. 
^As we have seen, the operator induced by bulk fields is localized only at y — ttR and thus the renor- 
malization prescription performed here is independent from the one of section B.2. 
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do not minimally couple to g_ , so that Sm'^^j^^ = 0. Summing over all colored fields, for 
periodic KK gluons (n > 0) we find 



»s 2 1 

m /„^ m 



A 



X 12. 



(B.5) 



We summarize in Table B.3 the different kind of contributions, summed over all the 
fields in the model. 





Amg 


i) bulk bosons 


27C{3) 


i) bulk fermions 


3 


ii) bulk bosons 


-f ln(2) 


ii) bulk fermions 





iii) boundary fermions 





Table 1: Summary of mass corrections in terms of 



For a cut-off scale A ~ (3-^4) /R, the mass splitting Arug turns out to be approximately 
equal to 

(B.6) 



Am„ = mf} - 2rrS^] ~ -1.4a^ — . 



C. Annihilation and coannihilation processes 

We collect in Table ^ all the matrix elements which are relevant for the computation of the 
DM relic density. Recall that the bulk fermions are in either the 81/3 or (iij^ of S'J7(3)^, 
where in the subscript we have denoted their U{1) charge under U{\)j^. After EWSB, 
they decompose as follows under SU{2)l x C/(l)y: 81/3 = 2i/g © I2/3 and 61/3 = 82/3 © 
2i/Q © l_i/3. In Table ^ we have denoted by Xi ^ and (f> respectively the SU{2)l singlet, 
doublet and triplet components of the lightest n = 1 KK mode of the 5D antiperiodic 
bulk fermions in both the 3 and the 6, with the understanding that for the 3 (j) (and 
the corresponding processes) are missing. These fields coincide with the states that we 
have collectively denoted by cl\ etc. in Fig. |l] and in the main text. The subscript 
a,b = 1,2 refers to the two distinct towers of KK mass eigenstates coming from the fermion 
pairs (^_,^'_). The SM fermions are denoted by / when we are considering both quarks 
and leptons and q for quarks only. We denoted by and the n = KK mode of 
the SU {2)l periodic triplets arising from the 5D bulk fermions ^_|!^, as in Fig. [l|. For each 
process, we also write the particle exchanged in the various {s, t, u) channels, whenever the 
flavour and gauge symmetries allow it. The channels mediated by should be considered 
only for the framework of Section 3.2. The fourth column 4p indicates when a four-point 
interaction vertex is present. 
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Diagrams 



Process 



t 



4p 



(1) aW 



~7ol — W 
9+ ,9+ 

(0) (2) 

9+ '9+ 

(0) 

9+' 



(0) 

9+ 



(0) (2) 

9+ ,9+ 

y+ '5+ 

(0) 

9+ 



(0) 

5+ 



9f\9f 



(0) 



£1) 



.(1) 



^(1) 

£1) 

.(1) 



(Xa,^a) 



XnXa - 
XaXa - 
XaXa 
XaXb 
XaXa 
XaXa,b 
- 

Ipa'ipa - 
i'a'tpa - 

ijjalpb — 
-Ipa-tpa - 
'>patpa,b 



4>a4>b 



(Pa(Pa,b - 
Xa,b'ipa,b 
^a,bi^a,b 
(t>a,b^a 



qnqn. 

qLQL 
(0) (0) 

9+9+ 
QMr 
• QlQl 

■ qmR 

■ bflf 

qLQL 
(0) (0) 

■ 9+9+ 

qLQL 

(0) (0) 

9+ 9+ 
bfbf 



9'y^ 



ly 
,(1) 



ly 
.(1) 



Ayi,a,b - 


(0) 

9 + m 

(0) 

5+ 9i 

(0),(0) 


Xa,6 




Xa,b 
^a,b 
4>a,b 




- 

9y9'y - 
9y9y - 


^ {qRm,qLqL) 
.bfbf 

(0) (0) 


(0) (2) 

5+ ,5+ 

(0) 

9+ 

(0) 


{Xa,'fpa) 
4>a 
9'y 


{Xa,tpa) 
4>a 

9^' 


X 








{Xa,i^a) 
4>a 


{Xa,i^a) 
4'a 




(1) 

9- Xa,b - 

9yi^a,b - 
9y<pa,b - 


(0) 

9+ QR 

(0) 

9+ qL 

(0)i,(0) 


Xa,6 
'^a,b 
4>a,b 


9'y 
9'y 
9y 


Xa,b 
'^a,b 
4>a,b 





Table 2: List of all the relevant (co-) annihilation processes. See text for details. 



State 


A^y 


9'^^ 


b'^y 






D.F. 


3 


24 




72 


48 



Table 3: Degrees of freedom for the states involved in coannihilation. 
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In Table | we list the degrees of freedom for the states relevant in the computation of 
the A^^^ relic abundance. For fermions we have D.F. = 2 x ANcUg, where is the color 
factor and Ug the number of states in the SU{3)w multiplet. The overall factor 2 takes 
into account the presence of two distinct towers for the antiperiodic fermions. In the case 
of gauge bosons one has simply D.F. = 3Ng, where Ng is the number of generators of the 
gauge group. 



D. Running of as 

The effective annihilation rate of our DM candidate is dominated by coannihilation and 
resonance effects with colored particles. As a result, the DM abundance is quite sensitive 
to the strong coupling constant ag'- ^dm oc a~^, where is here evaluated at the typical 
energy scale of the annihilation processes, namely 1/R. It is then important to evolve 
from, say, the scale of the Z boson mass mz up to 1/R taking into account the various 
threshold effects due to the many particle states in the model with a mass in that range. 
Due to the limited energy range of the RG flow and the presence of several states with 
a mass close to 1/R, we have computed the QCD /3-function at one-loop level in a mass 
dependent scheme (momentum subtraction) |^^, rather than in the more usual Minimal 
Subtraction (MS) or modified MS (MS) schemes. In this way, all the effects of threshold 
corrections are automatically taken into account, with no need of matching conditions. The 
drawback, of course, is that the resulting /3-function has a much more complicated form 
and the RG evolution has to be performed numerically. 

The computation is considerably simplified in a background field gauge (with ^ = 1 
gauge fixing parameter), where the /9-function is extracted by considering only the gluon 
2-point function. The contribution of a Dirac fermion and a gauge field with mass 
M in a representation r of SU{3)s can be written in a compact way as (5{gs, M/ fi) = 
gs{fi)Vl67T^[Pf{M/f,) + Pg{M/fi)], withio 

fi u^x'^(l-x)'^ 

P,(M/,) = STAr)lj^,^^^y (D.l) 

^ f^u'^x(l-x)(l+9x-6x'^) , 



and the convention that T[fund.) = 1/2. Eqs.( p.i| ) and ( p.2| ) give the correct contribu- 
tions for M ^ 0: I3f{g) 4/3r/(r) and f3g{g) Tg{r){-7/2 = -11/3 + 1/6), where 
— 11/3 and 1/6 are the contributions of the transverse and logitudinal (scalar) components 



of the gauge field. Eq. (D.l) is also in agreement with whereas we are not aware 
of the presence of an explicit expression for f3g in the literature. At one-loop level, one 



^"Notice that for unperturbed KK states (i.e. with vanishing bulk-to-boundary mixing terms), the (3- 
function contribution of the whole tower of KK states can be resummed, giving rise to /iii-dependent Coth 
or Tanh hyperbolic functions (for periodic or antiperiodic fields, respectively). One can then easily check 
that for _R ^ or _R ^ oo these functions interpolate between the usual 4D log-regime and the 5D linear- 
regime for Qs . One should however recall that in our model A ~ (3 4) / R so that for ^ ~ A we enter in a 
strong coupling regime where our perturbative expressions break down. 
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can consistently take the masses to be constant and /^-independent. Given the above 
relation, we have then numerically computed the RG evolution of as starting from its 
value as{Mz) = 0.117 up to 1/R, by including all particles in the model (for all states 
Tf = Tf{fund.) = 1/2 and Tg = Tgiadj.) = 3) with a mass up to 4/i?, to take also into 
account of residual threshold effects of more massive states. For /x < l/(2i?), ag decreases 
but then, in the range l/(2i?) < n < ^/R, due to the many colored fcrmions which become 
active, the sign of (3 changes and as starts to increase. The two effects tend to compensate 
each other, so that eventually as{l/R) is quite close to as{Mz)- For example, for 1/i? ~ 5 
TeV, as{l/R) c± 0.097 and as{l/R) ^ 0.092 for the model with one or two copies of SU{2>)s, 
respectively. Notice, just for comparison, that in a UED model with integer KK modes 
and 1/R ~ 2.5 TeV, we would have got a value significantly lower for ag'. as{2/R) ~ 0.075. 
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